Uniaxial dielectric mediums with hyperbolic dispersion 

relations 



Tom G. Mackay* 

School of Mathematics, University of Edinburgh, Edinburgh EH9 3JZ, UK 

Akhlesh Lakhtakia^ 
C ATM AS - - Computational & Theoretical Materials Sciences Group 

Department of Engineering Science and Mechanics 
Pennsylvania State University, University Park, PA 16802-6812, USA 

Ricardo A. Depine * 
Grupo de Electromagnetismo Aplicado, Departamento de Fisica, 
Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, 
Ciudad Universitaria, Pabellon I, 1428 Buenos Aires, Argentina 



Abstract 

The dispersion relations for conventional uniaxial dielectric mediums may be char- 
acterized as elliptical or elliptical-like, according to whether the medium is nondis- 
sipative or dissipative, respectively. However, under certain constitutive parameter 
regimes, the dispersion relations may be hyperbolic or hyperbolic-like. We investigate 
planewave propagation in a hyperbolic/hyperbolic-like uniaxial dielectric medium. For 
both evanescent and nonevanescent propagation, the phase velocity is found to be pos- 
itive with respect to the time-averaged Poynting vector. A conceptualization of a 
hyperbolic-like uniaxial medium as a homogenized composite medium is presented. 

Keywords: Hyperbolic dispersion relations, elliptical dispersion relations, Bruggeman ho- 
mogenization formalism 

1 Introduction 

As the materials sciences and technologies continue their rapid development, realistic pos- 
sibilities are emerging of realizing so-called metamaterials with novel and hitherto uncon- 
sidered optical/electromagnetic properties. A prime example is provided by the recently 
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discovered metamaterials which support planewave propagation with negative phase velocity 
(NPV), and thereby negative refraction. Until 2000, little attention had been paid to the 
phenomenon of negative refraction. Since 2000, there has been an explosion of interest in 
negative refraction [1, 2], following experimental reports of a metamaterial which supports 
negative refraction in the microwave regime [3]. 

Naturally-occurring uniaxial crystals have been extensively studied ever since the earliest 
days of the optical sciences. However, the electromagnetic properties of uniaxial mediums 
have recently been revisited by theoreticians in consideration of the prospects for NPV 
propagation in such mediums [4, 5, 6, 7]. A closely related issue concerns uniaxial dielectric- 
magnetic mediums with indefinite constitutive dyadics [8, 9]. 

The defining characteristic of a uniaxial dielectric medium is a distinguished axis of 
symmetry, known as the optic axis. Mathematically, the permittivity dyadic of a uniaxial 
dielectric medium may be expressed as 

| = el+ (e x - e) xx , (1) 

where a coordinate system has been selected in which the direction of the optic axis coincides 
with the direction of the unit vector x lying along the x axis, and J denotes the 3x3 identity 
dyadic. The real-valued parameter 

— for e x , e G R 

e 



7 



(2) 

Reie^} „ _ 
{ Reje} 



may be usefully employed to characterize planewave propagation in the medium specified 
by (1). The upper expression is appropriate to nondissipative mediums whereas the lower 
expression is appropriate to dissipative mediums. 

The electromagnetic/optical properties of uniaxial mediums with 7 > — this category 
includes naturally-occurring uniaxial crystals — have long been established. Comprehensive 
descriptions can be found in standard works [10, 11]. Uniaxial mediums with 7 < are much 
more exotic. Interest in these mediums stems from their potential applications in negatively 
refracting scenarios [8, 9] and in diffraction gratings [12], for example. 

Planewave propagation in a uniaxial medium is characterized in terms of a dispersion 
relation which is quadratic in terms of the corresponding wavevector components. The dis- 
persion relations for nondissipative mediums with 7 > have an elliptical representation, 
whereas a hyperbolic representation is associated with 7 < 0. In this communication we in- 
vestigate the planewave characteristics and conceptualization of uniaxial dielectric mediums 
with hyperbolic dispersion relations. 
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2 Planewave analysis 

The propagation of plane waves with field phasors 

K{r) = E_ exp (ik • r) 
Kix) — Ko exp (ik • r) 



(3) 



in the uniaxial dielectric medium specified by the permittivity dyadic (1) is investigated. 
The permittivity parameters are generally complex-valued; i.e., e, e x G C. The wavevector 
k is taken to be of the form 

k = ax + flz, (4) 

where a G M, (3 G C and z is the unit vector directed along the z axis. This form of k 
is appropriate to planar boundary value problems [11] and from the practical viewpoint of 
potential optical devices [12]. We note that the plane waves (3) are generally nonuniform. 
The source-free Maxwell curl postulates 



V x E(r) = iujB(r) 

V x H(r) = -iuoD(r) 



(5) 



yield the vector Helmholtz equation 

[(Vxl) • (Vxl)-^ 2 |] -E(r)=0, (6) 

with /i being the permeability of free space. Combining (3) with (6) yields the planewave 
dispersion relation 

(a 2 + p 2 — e/i cj 2 ) (a 2 e x + (3 2 e — e x efi cu 2 ) = . (7) 
In the following we consider the time-averaged Poynting vector 

= exp (-2 Im mz) R _ (g) 

Evanescent plane waves are characterized by Im{/9} > 0. The scenario characterized by 
Im{/5} < is not physically plausible for passive mediums and is therefore not considered 
here. 

2.1 Ordinary wave 

The ordinary wavevector 

k or = ax + (3 or z, (9) 
arises from the dispersion relation (7) with components satisfying 

a 2 + (3 2 or = cA/i . (10) 
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The vector Helmholtz equation (6) yields the eigenvector solution E_ = E y y, directed parallel 
to the unit vector y lying along the y axis, where the complex-valued magnitude E y is 
determined by the initial/boundary conditions. Consequently, the time-averaged Poynting 
vector reduces to 

£fe) = exp (-2 Im {&.};) (n) 

Since 

Re{A; or } • P(r) = z ) \ Ey \* [ a 2 + (Re {(3 or } f] > 0, (12) 

we say that ordinary plane waves have positive phase velocity (PPV) for all directions of 
propagation. 

Let us focus attention on a nondissipative medium (i.e., e, e x e R). From (10) we see that 
Im{/9 or } for (i) e > when u 2 efi < a 2 ; and (ii) e < 0. Thus, nonevanescent ordinary 
plane waves propagate in a nondissipative medium only when e > and —uj^/efi < a < 
Uy/e/j, . In geometric terms, the wavevector components have a circular representation in 
(a, P or ) space. 



2.2 Extraordinary wave 

The extraordinary wavevector 



hex = a £ + Pex 1 , 



arises from the dispersion relation (7) with components satisfying 



a e x + Pl x e 



uj e e x fi . 



(13) 



(14) 



In the case where f3 ex = the mathematical description of the extraordinary wave is iso- 
morphic to that for the ordinary wave. Therefore, we exclude this possibility from our 
consideration in this section. The eigenvector 



x — 



I E x , 



(15) 



arises as a solution to the vector Helmholtz equation (6); the complex-valued magnitude 
E x is determined by the initial/boundary conditions. The corresponding time-averaged 
Poynting vector is provided as 



P(r) 



exp(-21m{(3 ex }z) 



+ P*ex + " 




(16) 
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Hence, we find 



RefeJ .P(r) 



exp(-2Im{/3 ex }2;) 



a 



e(3 e 



+ Re 



(Re{/5 e:c }) 2 
ex/5* 



+ Re{ / 3 ez } Re 



e* 8* 



(17) 



We analytically explore the nondissipative scenario for nonevanescent and evanescent 
planewave propagation in Sections 2.3 and 2.4, respectively, whereas both the dissipative 
and the nondissipative scenarios are treated graphically in Section 2.5. 



2.3 Nonevanescent propagation 

By (14), the inequality 

u 2 e x fi - a 2 7 > (18) 

is satisfied for nonevanescent planewave propagation in a nondissipative medium, where 7 is 
defined in (2). Thus, lm{8 ex } = 0. We explore the cases 7 > and 7 < separately. 

(i) If 7 > then we require — cj^/e/i,, < a < ujy/efT in order to comply with (18). This 
implies that e > and e x > 0. In geometric terms, the wavevector components have 
an elliptical representation in (a, j3 ex ) space. 

(i) If 7 < then the inequality (18) reduces to cu 2 e/i < a 2 . Therefore, we see that 
nonevanescent propagation arises for (a) a > uj^/eji Q and a < —uj^fep^ when e > 0; 
and (b) —00 < a < 00 when e < 0. In geometric terms, the wavevector components 
have a hyperbolic representation in (a, f3 ex ) space. 

For lm.{3 ex } = and e X) e e K, we find that (17) reduces to 

RefeJ .P(r) = (19) 
Hence, nonevanescent plane waves have PPV regardless of the sign of 7 or e x . 



2.4 Evanescent propagation 

We turn to evanescent planewave propagation in a nondissipative medium as characterized 
by the inequality 

co 2 e x fx - a 2 7 < . (20) 

Hence, we have Re {f3 ex \ = 0. As in the previous subsection, we explore the cases 7 > and 
7 < separately. 
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(i) If 7 > then the situation mirrors that which we described earlier for hyperbolic 
nonevanescent propagation. That is, evanescent propagation arises for (a) a > uj^/efi 
and a < —uj^/efi when e > 0; and (b) — oo < a < oo when e < 0. In geometric terms, 
the wavevector components have a hyperbolic representation in (a, Im {{3 ex }) space. 

(ii) If 7 < then evanescent propagation arises provided that e > 0, e x < and —Uy/ejT < 
a < u^/efi . In geometric terms, the wavevector components have an elliptical repre- 
sentation in (a, Im {j3 ex }) space. 

For Re {f3 ex } = and e x , e e i, we find that (17) reduces to 

2 

RefeJ -P(r) = _;i exp(-2Im{/3 e ,}z) . (21) 

Hence, evanescent plane waves have PPV if (a) 7 < or (b) 7 > and e x > 0. However, 
negative phase velocity (NPV) propagation arises if 7 > and e x < 0. 



2.5 Illustrative examples 

Let us illustrate the geometric aspect of the dispersion relations with some representative 
numerical examples. 

First, suppose we consider the case 7 > with e = 2e and e x = 6e , where e is the 
free-space permittivity. In Figure 1 the real and imaginary parts of (3 ex are plotted against 
a. The elliptical nonevanescent nature of the dispersion relation is clear for —oo^/2t ^ < 
a < Uy/2e fj, , while the hyperbolic evanescent nature is apparent for a < — u; v / 2e /i„ and 
a > u^2e a fi . The elliptical/hyperbolic geometric interpretation breaks down when dissi- 
pative mediums are considered. However, the corresponding dispersion relations are geo- 
metrically reminiscent of their nondissipative counterparts. This can be observed in Fig- 
ure 2 in which the graphs corresponding to Figure 1 are displayed for e = (2 + i0.5) e and 
e x = (6 + ?0.75)e () . 

Second, we turn to the case 7 < with e = 2e and e x = — 6e . The real and imaginary 
parts of f3 ex are graphed against a in Figure 3. The graphs mirror those of Figure 1 but with 
nonevanescent and evanscent aspects interchanged; i.e., we observe hyperbolic nonevanescent 
characteristics for a < —uj^2e fi and a > uy/2e fi , and elliptical evanescent characteristics 
for — cj v / 2e () /i < a < u^2e fi . The corresponding graphs for e = (2 + i0.5) e„ and e x = 
(—6 + iO.75) e are presented in Figure 4. Notice that the shapes of the graphs in Figures 4 
and 2 are similar but not identical. 



3 Numerical conceptualization 

Although uniaxial dielectric mediums with 7 < do not occur in nature (to the best of the 
authors' knowledge), they can be conceptualized as metamaterials by means of homogeniza- 
tion. 
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For example, let us consider the homogenization of a composite comprising two compo- 
nent materials phases, labelled as a and b. Both component material phases are taken to 
be isotropic dielectric mediums: e a and e b denote the permittivity scalars of phases a and 
b, respectively. The component material phases are envisioned as random distributions of 
identically-oriented, spheroidal particles. The spheroidal shape — which is taken to be the 
same for all spheroids in component material phase a and b — is parameterized via the shape 
dyadic U_ = diag (U x , U, U). That is, we take the spheroid's principal axis to lie along the x 
axis. The spheroid's surface is prescribed by the vector 

L s (0A)=Vg'r(9, ( j ) ), (22) 

with r being the radial unit vector specified by the spherical polar coordinates 6 and 0. The 
linear dimensions of the spheroid, as determined by the parameter 77, are assumed to be 
small relative to the electromagnetic wavelength(s). 

The permittivity dyadic of the resulting homogenized composite medium (HCM) 

g HCM = diag (ef CM , e HCM , e HCM ) , (23) 

as estimated using the Bruggeman homogenization formalism, is provided implicitly via 

f a g? + fbi = g, (24) 

where f a and fb = ^ — fa denote the respective volume fractions of the material component 
phases a and b. The polarizability dyadics in (24) are defined as 

i = - f CM ) ' [L + iuR* {e e L - f CM ) ]~\ (t = a,b), (25) 
wherein the depolarization dyadic is given by the surface integral 

d = -^— [ d<p I de sine [ — J— = — 1 u- 1 -ff -tr 1 . (26) 

= iuAn J Y J yr -UT 1 -g HCM -UT 1 -rj= — = K J 

Closed-form expressions for the depolarization dyadic for uniaxial mediums are available in 
terms of hyperbolic functions [15]. However, we note that these exact results are not valid 
for nondissipative mediums with 7 < 0, and numerical evaluation of D has to be resorted 
to. 

The Jacobi iteration scheme 

f CM [p] = T {f CM [p - 1]} , (p= 1,2,3,...), (27) 

where the operator T is defined via 

T {g HCM } = {f a e a [L + iuR ' (e a L ~ f CM ) } 1 + f b e b [L + iuR • (e b L ~ f CM ) } 1 } 

.{/ a [L + iwR - (e a Z-f CM )] _1 + / 6 [L + iuR- (e^-f^)]- 1 }" 1 , 

(28) 
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may be employed to solve (24) for e HCM . Suitable initial values for the iterative scheme are 
provided by 

f M [0] = (// + /^)L (29) 

For further details on the Bruggeman homogenization formalism the reader is referred to 
[13, 14] and to references therein. 

Let us consider the homogenization scenario wherein material component phase a is taken 
to be iron at 670 nm free-space wavelength. Correspondingly, we take e a = (—4.34 + i20.5) e . 
The material component phase b is assumed to be free space; i.e., e b = e . The shape of 
the component spheroids is specified by U x /U = 12. The Bruggeman estimates of the HCM 
permittivity parameters e HCM and e^ CM are plotted as functions of volume fraction f a in 
Figure 5. At intermediate values of f a we see that 7 < for a substantial range of f a values. 

Extensive accounts of similar numerical homogenizations, based on the Bruggeman for- 
malism and more general approaches, can be found elsewhere [13, 16]. 

4 Concluding remarks 

The dispersion relations for uniaxial dielectric mediums have been characterized with respect 
to the parameter 7 (2). For 7 < 0, the dispersion relations are hyperbolic for nondissipative 
mediums and hyperbolic-like for dissipative mediums. Similarly, the dispersion relations 
are elliptical for nondissipative mediums and elliptical-like for dissipative mediums with 
7 > 0. Through the homogenization of isotropic component material phases based on 
spheroidal topology, we demonstrate that metamaterials with 7 < may be straightforwardly 
conceptualized. Thus, a practical means of achieving the exotic electromagnetic properties 
associated with hyperbolic and hyperbolic-like uniaxial mediums is presented. 
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Figure 1: A plot of the real (solid curve) and imaginary (dashed curve) parts of (5 ex against 
a for e x = 6e and e = 2e . The values of a and f3 ex are normalized with respect to Uy/e fj, . 




Figure 2: As figure 1 but for e x = (6 + i0.75) e and e = (2 + i0.5) e . 
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Figure 3: A plot of the real (solid curve) and imaginary (dashed curve) parts of f3 ex against 
a for e x = — 6e and e = 2e . The values of a and f3 ex are normalized with respect to Uy/e fi . 
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Figure 4: As figure 3 but for = (—6 + i0.75) e and e = (2 + i0.5) e . 
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Figure 5: The real (above) and imaginary (below) parts of e HCM and e^ CM plotted against 
volume fraction f a . The permittivity values are normalized with respect to e . Component 
phase values: e a = (—4.34 + ?20.5) e and e b = e ; spheroidal shape parameters: U x /U = 12. 
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